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to  attract  (Cnntinu*  on  rmvrao  atdm  It  nocoooary-  and  tdmnltly  by  block  nufnh^,} 

Efficient  compressiion  of  low  mass-density  payloads  by  the  implosion  of  higher  ma's-density 
liquid  cylinders  or  liners,  as  in  the  NRL  LINUS  concept  for  controlled  thermonuclear  fusion, 
requires  rotation  of  the  liner  material  to  avoid  Rayleigh -Taylor  instabilities  at  the  liner-payload 
interface.  Experimentally,  such  implosions  have  been  demonstrated  with  liners  formed  within 
rotating  implosion  chambers.  The  pre.sent  work  uses  a scale-model  experimental  apparatus  to 
investigate  the  possibility  of  creating  liner  rotation  by  tangential  injection  of  the  iiquid  liner 
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^20.  Abstract  (Continued) 

material.  Different  modes  of  behavior  are  obtained  depending  on  the  fluid  exhaust  procedures. 
Right<irrular,  cylindrical  free  surfaces  are  achieved  with  axial  exhaust  of  fluid  at  radii  interior 
to  the  injection  nozzles,  for  which  the  liner  exhibits  a combination  of  solid-body  and  free  vortex 
flows  in  differtMU  regions.  Measurements  allow  estimates  of  power  losses  to  viscous  shear,  turbu- 
lence, etc,  A simple  model  based  on  open-channel  flow  is  then  derived,  which  is  in  good  agreement 
with  experiment,  and  is  used  to  extrapolate  results  to  the  scale  of  a possible  LINUS  fusion  reactor. 


CONTENTS 


I.  INTRODUCTION 


II.  DESCRIPTION  OF  THE  APPARATUS 


III.  FLOW  DESCRIPTION 


IV.  experimi:nt 


V.  SWIRL  CHAMBER  FU1W  MDDEI 


POWER  LOSS  CALCULATIONS 


CONCLUS IONS 


VIII.  REFERENCES 


AcQosslonFor  y ] 

NTIS  GKAAI 

DDC  TAB 

1 

UncU'.nour.o?d 

j Justification 

! Pr  ■ ■ 

i 

..1 

C.KNEKATION  OF  A ROTATINCJ  LIQUID  LINER  BY  TANGENTIAL  INJECTION 


I.  INTRODUCTION 

Plasma  heating  by  imploding  liquid  liners  is  being  investigated, 
at  the  Naval  Research  Laboratory,  as  an  attractive  concept  for  a 
thermonuclear  fusion  reactor.  In  the  LINUS  concept  (1),  a stable 
D-T  plasnui  is  created  in  the  cavity  formed  by  rotating  a liquid 
metal  cylinder  or  liner  inside  a partially  filled  pressure  vessel. 

A cylindrical  free  surface  is  thus  formed  at  a radius  surround- 

ing a plasnw  and  magnetic  field.  The  liner  rotation  must  be 
sufficient  to  prevent  Raylelgh-Taylor  instability  of  the  inner 
free  surface  of  the  liner  during  final  compression  of  the  plasma 
and  magnetic  field  (2,3).  The  liner  is  cojtstrained  on  its  outer 
surface  by  a free  piston,  which  can  be  rapidly  accelerated  by  a high 
pressure  gas  driver.  The  resulting  free  piston  displacement  produces 
an  inward  radial  motion  rj(t)  of  the  free  surface,  which  at  high 
magnetic  Reynolds  number  produces  an  adiabatic  compression  of  the 
magnetic  field  and  plasma,  raising  the  plasma  temperature  above 
thermonuclear  ignition.  The  resultivig  energetic  neutrons  are 
captured  in  the  liner  to  raise  its  temperature,  react  with  Li  to 
breed  tritium,  and  protect  the  reactor  structure  from  radiation 
damage.  This  process  in  a reactor  is  typically  repeated  at  a rate 
of  a few  Hz . 

A number  of  experiments  and  theoretical  studies  of  the  LINUS 
concept  have  been  performed  (1,  3-7).  In  all  laboratory  experi- 
ments to  date,  the  initial  liner  configuration  has  been  created 

by  rotating  the  entire  pressure  vessel  (Implosion  chamber)  , 

NoU*;  ManuiirripC  «uhmitl<'<l  April  19,  1979. 
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so  that  iho  is  in  solid-body  rotation  when  the  implosion 

begins.  11  this  approach  is  carried  to  the  reactor  stage,  rotation 
ot  the  implosion  chamber  implies  that  a number  of  difficult  engineering 
problems  must  be  adviressed,  such  as  design  of  support  bearings  and 
bearing  lubricant,  and  design  of  a rotating  vacuum  seal  and  a high 
pressure  seal  for  the  pisti'n  gas  drive. 

Simple  heat  bal.ince  calculations  on  reactor-sized  systems 
producing  power  at  a level  of  several  hundred  MW(t)  show  that  a 

3 

liner  volumetric  flow  rale  ol  several  m /sec  is  required,  and  that 
this  must  be  supplied  by  external  pimips,  directly  coupled  to  a set 
of  exit  turbines  to  recover  tlie  considerable  discharge  kinetic  energy 
ol  the  liquid  metal.  II  these  same  pumps  covild  be  used  to  provide 
liner  rotation,  then  the  above  mentioned  engineering  com|>l Icat ions 
associated  with  a rotating  implosion  chamber  could  be  avoided.  A 
sm.ill  sc. lit'  exiH'riment  is  described  here  in  which  a rotating  liner  is 
treated  by  tangent i;i  I injection  of  liquid,  and  measurements  are 
m.ide  to  determine  the  power  consumed.  A simple  channel-flow  model 
is  Invoked  to  explain  the  power  loss  by  turbulent  shear  at  tlie 
chamlier  walls,  .ind  to  estim.ite  the  power  consumed  in  a reactor 
appl icat ion. 

1 1 . lIKSt^irriON  THK  ArPAKATUS 

In  Fig.  I,  a schenvit  Ic  is  shown  of  the  cylindrical  vessel  or 
'swirl  ch.imher'  with  tangent  i ;il-lnject  ion  nozzles  at  the  outer  radius 
.utd  exit  ports  mounted  in  the  side  wall.  The  swirl  chamber  has  a 
r.idlus  i>t  12.5  cm  and  a length  of  b.8  cm.  Four  linear  tangential- 


lujoction  nozzles  are  provided,  eaot\  with  ii  rectangular  0.68  cm 


throat  of  dimensions  6.8  cm  x 1.0  mm.  The  nozzles  are  fed  from  a 


manifold  encircling  the  chamber  through  four  1.76  cm”"  lu’ses, 

.o 


which  enter  radially  at  the  outer  wall.  The  flow  turns  90  into 
a 6.8  cm  X 3 mm  settling  chamber  upstream  of  the  nozzle,  and  then 
discharges  along  the  outer  wall. 

A variety  of  exhaust  ports  are  provided  In  the  ch;unher.  A 

4 cm  diameter  port  at  the  center  exposes  the  free  surface  to 

atmospheric  pressure  at  all  times.  Four  sets  of  twin  discharge 
2 

ports  of  1.33  cm*" , are  located  midway  between  the  Injection  nozzles 
in  the  outer  wall,  (Fig.  1)  permitting  radial  discharge  flow.  Those 
ports  connect  through  eight  1.5  cm  diameter  hoses  to  a valved 
manifold,  permitting  the  fl(.>w  rate  to  be  controlled.  Finally,  a 
ring  of  8 ports,  1.3  cm  in  diameter  and  placed  at  10  cm  radius,  are 
located  in  one  end  wall  of  the  chamber.  These  ports  are  connected 
through  1.5  cm  diameter  hoses  to  a common  valved  m.anifold,  allowing 
the  flow  rate  to  be  controlled. 

The  end  wall  opposite  these  exhaust  ports  is  made  of  Plexiglas 
so  that  photographs  of  the  swirling  flow  can  be  taken.  A static 
pressure  gauge  in  the  outer  wall  reads  the  rotational  head  of  the 
spinning  fluid,  and  a pressure  gauge  on  inlet  manifold  gives  the 
source  pressure.  The  rotational  speed  of  the  free  surface  was 
measured  by  seeding  the  flow  with  3 mm  diameter  styrofoam  balls 
which  floated  on  the  surface  and  were  then  illuminated  with  a 
variable-frequency  stroboscopic  light  source. 
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III.  FLOW  DESCRIPTION 


The  initial  concept  was  to  inject  fluid  tangentially  and  extract 
it  at  the  outer  wall,  so  that  the  liquid  flow  interior  to  a thin 
layer  at  the  outer  wall  would  approximate  solid  body  rotation. 
However,  it  was  not  possible  by  this  method  to  produce  a circular 
free  surface.  Instead,  the  inner  surface  of  the  fluid  formed  a 
twisted  cylinder  of  quasi-elliptical  cross-section  whose  major  axis 
varied  in  orientation  by  90°  from  one  endwall  to  the  other.  This 
non-circular  condition  is  remedied  by  allowing  liquid  to  discharge 
axially  through  the  chamber  end  wall,  until  approximately  50%  is 
discharging  radially  and  50%  axially.  This  produces  a turbulent 
but  steady  free  surface  of  circular  cross  section.  As  the 
percentage  of  axial  discharge  is  increased  and  the  radial  discharge 
is  correspondingly  decreased,  the  free  surface  remains  circular 
and  steady.  The  chamber  can  be  operated  with  50  to  100%  of  the  flow 
discharging  axially.  The  case  of  100%  axial  discharge  provides 
the  simplest  flow  pattern,  so  that  the  radial  discharge  at  the 
outer  surface  was  reduced  to  zero  for  the  experiments. 

At  the  maximum  flow  rate  of  3.3  1/sec  the  diameter  of  the  free 
surface  is  14.8  cm,  which  is  smaller  than  the  location  of  the  axial 
exhaust  ports.  This  means  that  the  fluid  between  the  free  surface 
and  these  axial  ports  should  nearly  be  in  solid  body  rotation. 

The  flow  pattern  between  the  axial  exhaust  ports  and  the  outer  wall 
is  more  complicated,  and  can  be  divided  into  two  regions.  In  the 
inner  region,  flow  is  rotational  and  possesses  an  inward  radial 
velocity,  similar  to  a free  vortex.  In  the  outer  region  the  flow 


h 


is  azlmuthally  periodic,  having  a high  tangential  velocity  at  the 
nozzle  exit,  (and  zero  velocity  at  the  radial  exhaust  ports, 
if  these  are  in  use). 

IV.  EXPERIMENT 

With  a circular  free  surface  established  in  the  chamber,  the 
following  quantities  were  measured  (Table  I): 

Table  I 


Total  Injected  flow  rate, 

Supply  pressure: 

Outer  wall  static  pressure,  p 

w 

Free  surface  pressure,  p^: 

Free  surface  radius,  r ^ : 

Free  surface  rotational  speed,  : 


3.3  1/sec 

2.43  atm  (gauge) 
0.32  atm  (gauge) 
0.00  atm  (gauge) 

7.4  cm 

63  radians/sec 


From  this  data  the  following  quantities  are  calculated  (Table  II). 


Table  II 

Nozzle  throat  injection  velocity, 

Exhaust  flow  velocity,  V : 

c. 

Free  surface  rotational  velocltv,  : 

”l 

2 

Fluid  power  injected,  SmV  jjjji 

, 2 

Fluid  power  exhausted,  : 

E 


11.3  m/sec 
3.5  m/sec 
4.7  m/sec 
212  watts 
20  watts 


V.  SWIRL  CHAMBER  FLOW  MODEL 

In  order  to  calculate  the  relative  Importance  of  various  power 
dissipation  mechanisms  in  the  swirl  chamber,  a velocity  distribution 
consistent  with  the  experiment  must  be  estimated.  The  flow  is  assumed 


S 


to  bi“  inJopoiuloiU  ot  axial  position,  aiul  is  d iv ivlid  radially  into 


tlirof  ro^ions  as  disoussod  abovo.  In  tin-  innor  ri'^ion,  insido  the 

radius,  I , ot  the  axial  exhaust  duets,  the  lailial  veloeitv  is  aero, 
r 

This  t low  is  taken  to  be  in  si>lid  bovly  rotation,  with  the  inner  free 
surtaee  rotatln^t  at  4.7  m/see.  Kor  radii  larger  than  the  exh.iust  ports 
the  tlow  divides  into  two  ret;ii>ns.  The  injected  fluid  torms  <i  radially 
thin  i»'t  at  the  outer  wall,  anil  mixes  turbulent  Iv  with  an  inner, 
.i/imuthally  uniform  flow  which  beh.ives  as  .1  t re*-  vortex. 

The  inner  solid  body  a::imutha!  velocity  is 
V - V^r/r^  r ^ r^,. 

where  Vj  is  the  me.isured  .i/imuthal  free  surface  velocity.  Conserving 
an>;ular  mi'mentuni,  Vr  , in  the  tree  Vv'rtex  flow  };ives 


V^r-^/r^r 


'"V. 


The  free  vortex  flow  must  now  be  matched  to  the  lnJectioi\  region. 

This  is  done  in  a consistent  manner  by  assuming  that  the  Jet  flow 

from  the  tangential  injection  nozzles  expands  radially  inward,  such 

that  the  mean  Jet  velocity  equals  the  free  vortex  velocity.  Thus, 

if  the  jet  thickness  Is  C)r  and  the  nozzle  width  is  ar  , the  radius 

r*  - r -Or  at  which  the  let  and  free  vortex  flows  match  is  defined 
w 


bv; 


V* 


IN.I 


\ V.r,,*" 

V*-*  ‘i*-* 


J 


Sv'lving  for  r*: 


and  for  6r: 


Or  ■ r -r*  - r (1- 

V w 


— ) 


For  4^r  - 0.1  cm,  6r  - 0.225  cm,  so  r*  - 12.275  cm  and  V*  ■ 5.17  m/sec. 

The  mean  velocity  on  the  outer  wall  will  lie  between  V*  and 
ViN).  Transit  ion  from  to  V*  will  occur  over  a Jet  length 

estimated  by  assuming  that  the  Jet  Inner  boundary  is  parallel  to 
streamlines  of  the  free  vortex  flow.  The  radial  velocity  U*  is 


2nr*h 


" 6.26  cm/sec 


The  streamline  slope  is  S • ll*/V*  » 0.0121,  from  which  the  Jet 
length  is  L - (6r-*iR)/S  - 10.3  cm,  corresponding  to  an  angular  tran- 
sition distance  of  47.4  degrees.  For  a linear  variation  of  velocity 

during  transition,  the  mean  speed  at  the  wall  is  V(r  ) - 6.78  m/sec. 

w 

The  complete  radial  distribution  of  V can  be  modelled  as  shown  in 
Fig.  2 . The  velocity  distribution  of  Fig.  2 can  then  be  integrated 

Co  find  Che  azlmuChal  flow: 


■/ 


hV(r)dr 


wtiich  >;lvea  “ 19.4  1/sec. 

V 1 • _U’SJj  CALCULATIONS 

The  injected  fluid  power  (Table  II)  is  212  watts.  This  power  is 
distributed  to  exhaust  flow  power  (20  watts),  pressure  drop  at  the 
entrance  to  the  exhaust  duct,  turbulent  mixing  of  the  jet,  and  turbu- 
lent shear  at  the  wall. 

0 

The  exiiaust  duct  pressure  drop  is  approximately  equal  to 
for  sharp-edged  entrance,  giving  -‘ii,,  » 0.124  atm.  The  power 


f 


distribution  (Fig.  2)  we  get  (V  = 5.63  m/sec.  The  Reynolds 
~~2  h 

number  based  on  (V  ) and  m is  Re  = 11,500,  which  for  a relative 
roughness  of  .002  gives  f = .033  (9).  Thus  the  hydraulic  slope  is 
S - .669.  The  equivalent  pressure  drop  over  the  channel  length  ^ 
is  iipg  = Pg-^S,  where  ^ is  taken  as  the  circumference  at  the  mean 
channel  radius.  The  calculated  value  is  ^p  = .0412  atm.  The  power 
dissipated  in  wall  shear  is  thus  calculated  to  be  P.  = Qq^Pq  = 80  watts. 
The  calculated  power  loss  mechanisms  are  shown  in  Table  III: 


Table  III 


Calculated  Power  Losses  in  the  Swirl  Chamber 


Exhaust  Flow  Power 


20  watts 


Exhaust  Duct  Pressure 
Drop 


41  watts 


Turbulent  Jet  Mixing 


62  watts 


Wall  Shear  Loss 


80  watts 


Calculated  Loss: 


'03  watts 


Measured  P, 


212  watts 


The  calculated  loss  accounts  for  the  measured  power  input 
to  within  4%,  which  must  be  regarded  as  good  agreement  considering 
the  approximate  character  of  the  model  used.  For  example,  the  open 
channel  flow  model  underestimates  the  dissipation  at  the  outer 
radius,  where  the  velocity  is  considerably  higher  than  the  mean 
value  of  V used  in  the  calculations. 

Extrapolation  to  LINUS  Fusion  Reactors 


In  the  case  of  a LINUS  fusion  reactor  system  (8) , the  initial 
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Innor  surlaot!  railiiis  and  tannoiuial  speod  would  bo  about  one  moler  and 
20  m/soc,  respect Ivoly,  usln^  a liner  of  lead-lltbium  alloy  at  about 
700-800*  K,  (p  •“  10*  kj;/m  ^) . Witli  fixed  flow  neometry  (and  the  same 
relative  rouj’lmess  factor  of  0.002),  the  power  losses  in  a reactor- 
size  system  can  be  scaled  from  the  present  swirl  chamber  experiment 
as  sliown  in  Table  IV: 


Table  IV 


ScaKni  Kstlmates  ol  Power  bosses  In  blNUS  Reactor 


bxbaust  Ouct  Pressure  Drop 


().9  MW 


Turbulent  -let  Mlxiitp 


10.5  MW 


W.ill  Shear  boss 


9.4  MW 


Total  boss 


2 7 


Vor  a blNflS  t us  Ion  reactor  scenario  (as  discusseil  in  Ref.  8),  with 


a tliermonuc leaf  ener>;y  nain  relative  to  total  enerpy  of  1.5,  and  a 
peak  operating,  map.nctic  lleld  ot  0.8  MO,  a module  based  on  a simple 
scale-uii  ot  the  swiri  chamber  wouUl  provlile  about  216  M.l  ot  thermal 
ener>;y  or  about  72  Md  of  electrical  ener>;y  eacit  shot.  In  order  for 
the  flow  power  losses  to  be  U'ss  than  ten  percent  of  the  electrical 
output  of  tlu‘  reactor,  the  Implosion  repetition  rate  nuist  exceed 
3.7  Hz.  The  net  electrical  power  per  module  would  tlien  be  about 
240  MW.  Under  these  same  ci'nditions,  the  volume  flow  rate  of  liner 
material  ft>r  a mean  t emper.-it  ui't'  excursion  of  100'\'  Is  4.5  mVsec, 

3 

compared  to  a scaled-np  injection  Mow  rate  cvf  2.7  m /sec  to  proviile 
liner  rotation.  Thus,  power  losses  and  turbulence  levels  associated 


with  liner  circulation  to  limit  the  temperatxire  excursion  of  the 


“w  < 


1 

I 


liner  material  can  be  comparable  to  those  incurred  in  using  tangential  i 

injection  to  create  liner  rotation.  Separate  study  (10)  has  indicated 
that  fluid-bearing  losses  involved  in  rotating  the  entire  implosion 
chamber  are  comparable  to  those  estimated  for  the  tangential  injection 

I 

approach.  Thus,  considerable  gain  in  simplicity  of  LINUS  reactor  [ 

» 

operation  may  be  possible  without  incurring  significantly  worse  f 

i 

performance  by  use  of  tangential  fluid  Injection  to  form  the  rotating  f 

liquid  liner. 

I 

vn.  CONCLUSIONS  ; 

Experimental  tests  with  the  swirl  chamber  demonstrate  the  formation  , 

t 

of  rotating  liquid  liners  by  tangential  injection  of  fluid.  Cylindrical  ' 

inner  surfaces  with  satisfactory  circularity  and  axial  uniformity  are 
achieved  with  tangential  injection  at  the  outer  radius  of  the  chamber 

I ; 

and  substantial  ()50%)  extraction  of  the  fluid  mass  flow  at  inner  radii  ii 

through  ports  in  the  sidewalls.  A model  of  the  azimuthal  and  radial  j 

fluid  velocity  distributions  provides  a reasonable  basis  for  estimates  !■ 

of  power  losses  in  terms  of  an  open-channel  shear  flow  and  turbulent  p 

entrance  and  exit  conditions.  The  calculated  power  loss  is  in  good 

agreement  with  experimental  measurement  lending  credence  to  the  flow  ^ 

1, 

model  and  allowing  extrapolation  of  calculations  to  larger  devices.  I 
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